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Abstract: 

The nonequilibrium phase transition has been studied by Monte Carlo simulation 
in a ferromagnetically interacting (nearest neighbour) kinetic Ising model in presence 
of a sinusoidally oscillating magnetic field. The ('specific-heat') temperature derivative 
of energies (averaged over a full cycle of the oscillating field) diverge near the dynamic 
transition point. 



*E-mail:muktish@physics.iisc.ernet.in 



1 



I. Introduction 

The dynamics of Ising system in presence of an oscillating magnetic field has been 
studied extensively [1,3] in the last few years. The dynamic hysteretic response [1,3] and 
the nonequilibrium dynamic phase transition [2,3] became an interesting and important 
part of these studies. The nonequilibrium dynamic phase transition in the kinetic Ising 
model in presence of a sinusoidally oscillating magnetic field has been studied first by 
Tome and Oliviera [2]. They have solved the dynamic mean field equation of motion 
for the average magnetisation of the kinetic Ising model in presence of a sinusoidally 
oscillating magnetic field. Defining the dynamic order parameter as the time averaged 
magnetisation over a full cycle of the oscillating field, they showed that in the field 
amplitude and temperature plane there exists a phase boundary separating dynamic 
ordered and disordered phase. Precisely, below the phase boundary the value of the order 
parameter is nonzero and it gets zero value above the bounday. There exists a tricritical 
point on the phase boundary line separating the nature (discontinuous/continuous) of 
the dynamic transition. 

Later, Acharyya and Chakrabarti [3] studied this nonequilibrium phase transition 
in the kinetic Ising model by extensive Monte Carlo simulations. They have also studied 
[3] the temperature variations of the AC susceptibility components. Their important 
observation was that the imaginary (real) part of AC susceptibility gives a peak (dip) 
at the dynamic transition point, indicating the thermodynamic nature of this type of 
nonequlibrium dynamic phase transition. 

Here, in this paper the time averaged (over a full cycle of the oscillating field) coop- 
erative energy and total energy of the system have been calculated by Monte Carlo sim- 
ulation and studied the temperature variation of the temperature derivatives ('specific- 
heat') of those energies. These 'specific- heat' like quantities are observed to diverge very 
near the dynamic transition point (where the dynamic order parameter vanishes). The 
paper is organised as follows: in section II the model and the simulation scheme are 
discussed very briefly, the results are reported in section III and the paper ends with 
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some concluding remarks in section IV. 



II. Model and Simulation 

A ferromagnetically interacting (nearest neighbour) Ising system in presence of an 
oscillating magnetic field can be represented by the Hamiltonian, 

<ij> i 

Here, sf(= ±1) represents Ising spin variable and, h(t) = hocos(ujt) is the externally 
applied oscillating magnetic field, ho and u are the amplitude and frequency of the 
field. For simplicity all Jij(> 0)'s are taken here equal to unity and boundary condition 
is periodic. 

Here, a square lattice of linear size L(= 100) has been considered. At any finite 
temperature T and for a fixed frequency (a;) and amplitude (ho) of the field, the dy- 
namics of this system has been studied here by Monte Carlo simulation using Glauber 
single spin-flip dynamics. Each lattice site is updated here sequentially and one such 
full scan over the lattice is defined as the time unit here. 

The following quantities are calculated, for a fixed values of T, uj and ho'- 

(i) The dynamic order parameter, Q = (u/2-k) § m(t)dt, where m(t) is the instan- 
teneous magnetisation. This is essentially the time (t) averaged magnetisation over a 
full cycle of the oscillating magnetic field. 

(ii) The time averaged (over a full cycle) cooperative energy of the system 



Ecoop 



-(u/2nL 2 )j(^sls^jdt. 



(iii) The time averaged (over a full cycle) total energy of the system 

E tot = -(u/2kL 2 ) I I + h (t) £ si J dt. 

J \<ij> i J 

Each data point has been obtained by averaging over 50 different random Monte 

Carlo samples. 
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III. Results 

The temperature derivatives of E coop and E tot can be called 'specific-heat' for this 
nonequilibrium problem. The temperature variations of Q , C coop = 5E coop /5T and 
Ctot = 8Etot/8T have been studied. For a fixed values of the field amplitude ho, the 
dynamic order parameter Q shows a continuous (depending upon the value of ho) phase 
transtion at temperature Td(ho) [3]. The values of ho (=0.4 and 0.8) are chosen here in 
such a way that Q always undergoes a continuous transition. The temperature variations 
of Q, C coop and C to t have been shown in Fig. 1. In this case the frequency (uj) of the 
field is kept fixed (a; = 0.0628). From the figure it is clear that the 'specific-heat' s C coop 
and Ctot diverge near the dynamic phase transition point. Ctot shows more prominent 
divergence than that showed by C coop - Some of the results have been checked for larger 
sizes (L = 200) and no significant deviation was observed. The CPU time used here, to 
generate the whole set of data in IBM workstation, is approximately 16 hours. 

IV. Concluding remarks 

In conclusion, the nonequilibrium phase transition in the kinetic Ising model has 
been studied by Monte Carlo simulation. The 'specific- heat' shows divergence near the 
dynamic transition point. The divergence of 'specific- heat' near the dynamic transition 
point is a distinct signal of phase transition and an indication of the thermodynamic 
nature of this class of nonequilibrium phase transition. The indication of the thermody- 
namic nature of this type of nonequilibrium transition was first observed by Acharyya 
and Chakrabarti [3] in their study of the temperature variation of AC susceptibility. 
This observation of 'specific- heat' divergence near the transition point re-confirms, by 
another independent way, the thermodynamic nature of dynamic phase transition. A lot 
of computational effort is required to establish precisely the divergence of 'specific- heat' 
near the transition point. The work is in progress in this front and the details will be 
reported elsewhere. 

It would be quite interesting to know how the dynamic order parameter goes to 
zero and specific heat diverges as one approaches the transition point. The universality 
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class of this type of nonequilibrium transition, is not yet known. 
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Figure Caption 

Fig.l. Temperature variation of Q (solid line), C coop (circle) and C tot (triangle) for 
fixed to (= 0.0628) and two different values of the field amplitude (ho): (I) h = 0.4 and 
(II) ho = 0.8. L = 100 here. 
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